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LIFT AND MOMENT EQUATIONS FOR OSCILLATING AIRFOILS 
IN AN INFINITE UIKTAGGERED CASCADE 
By Alexander Mendelson and Robert ¥. Carroll 


SUMMARY . 

Exact equations are derived for the oscillatory aerodynamic forces 
acting in an unsta^gered cascade of airfoils fluttering in potential 
flow. Aerodynamic coefficients similar to those of the isolated airfoil 
are obtained as functions of the cascade geometry and the phasing bet'vreen 
successive blades j the phasings considered are zero, 90°, and 180°. It 
is shown that 90° is a special case of 180° phasing. These- aerodynamic 
coefficients are plotted for the special case when all the airfoils are 
vibrating in bending in phase (360° phasing) . It is shown that the ef- 
fect of cascading for this case is to reduce greatly the aerodynamic 
dan^jing. 


INTROIODTION 

The flutter of airfoils in a cascade has until recently been pri- 
marily of academic Interest. However, the widespread use of con 5 >ressors 
and turbines in current aircraft power plants has given the problem sig- 
nificance. Conpressor blades, in particular, are susceptible to vibra- 
tions, and some of these vibrations have been attributed to flutter. 

The problem of the flutter of a con 5 )ressor or turbine blade differs 
from that of an Isolated airplane wing in at least two ways. It is ne- 
cessary to consider, first, the effect of centrifugal force; and, second, 
the effect of cascading. The effect of centrifugal force can be tahen 
into account with sufficient accxiracy by applying the appropriate centri- 
fugal force correction to the fundamental bending frequency (ref. l). 

The effect of cascading is much more difficult to evaluate, however, 
since it is necessary to take into aecount the interference effect be- 
tween the blades, which obviously depends on the cascade geometry. Two 
new geometric variables must therefore be Introduced; namely, the spaning 
between blades and the stagger angle. In addition, since flutter is a 
time-varying phenomenon, another parameter must be introduced to take 
accoTint of the phasing between the motions of the different blades of 
the cascade . 
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The problem of flutter of coiipressor blades is first simplified by 
assuming an infinite cascade of aixfoils. A first step in a flutter 
analysis of such a cascade of airfoils is to determine the oscillatory 
aerodynamic forces and moments acting on the cascade. Most recorded 
cases of compressor blade fliitter Indicate the occurrence of such flut- 
ter at hi^ aerodynamic loading, -where the blade stalls and flow separa- 
tion occurs, and it woiild be desirable to solve the problem for this 
case . However, no general methods _ for calculating aerodynamic forces in 
nonpotential flow are available. It is therefore necessary, as a first 
approach, to consider the case in potential flow at low angles of attack. 
The effect of flow separation at s-ba.11 can then be taken into account 
separately, for instance, by the introduction of aerodynamic time lags 
(refs. 2 and. 3) or other mechanisms -which may prove -useful. The object 
of this paper is to present solutions for several special cases of the 
oscillatory aerodyna mi c forces and mom^ts acting on an infinite cascade 
of airfoils in potential flow. 

The first derivation of oscillatory forces in a cascade was made in 
reference 4. The effect of the wind-tunnel -walls on a fluttering isolated 
airfoil -was determined approximately. This is equi-valent to the specieil 
case of an infinite cascade -without stagger, -vrLth adjacent blades being 
180° out of phase. The integral equation for the problem -was set up and 
an approximate solution obtained by replacing the kernel with a simple 
polynomial. The results are not applicable for spacing-to -chord ratios 
of less than 1. The same problem was solved rigorously in reference 5 
and in reference 6 by different methods. The results are obtained in the 
form of do-ubly infinite series of Jacobian elliptic functions. 

Another special case -was treated in reference 7. A cascade -with 
stagger -with al l the blades -vibrating in phase is considered. Hie method 
is similar to that xised in reference 8 for the isolated airfoil. The 
form in -which the final results are presented cannot be easily used for 
numericeLl calculations. More recently, the general integral equation for 
a cascade with stagger and prescribed phasing -was set iq) by Sisto (ref. 

9) . Approximate numerical solutions -were then obtained for se-veral cases 
of zero stagger. The method -used is similar to that of reference 4. 

The present paper attenpts to fill some of the gaps left by the pre- 
-vlo-usly mentioned investigations. A solution is obtained for the case of 
zero stagger -with the blades either 180° out of phase (considered in ref. 
6), 90° out of phase, or in phase. The last -two cases can be handled 
approximately by the method of reference 9. The present paper presents 
an exact solution. Furthermore, the results are obtained in a form which 
allows a solution in a sinqjle s-tiai^tforward manner on a desk calculator, 
the final results being presented in the form of aerodynamic coefficients 
-which can be calculated in a step-wise manner by "the use of a set of recur- 
rence formulas. 
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AERODYMMIC LIFT AND MOMENT 

The analysis of the oscillatoiy forces and moments acting on a cas- 
cade of airfoils is made using the classical methods. In the beginning, 
the cascade is assumed to he of infinite extent with arbitrary spacing 
and stagger, as shown in figure 1. The problem is later specialized to 
the case of zero stagger. The airfoils are assumed to he thin and per- 
forming small oscillations in a potential, inccai 5 )ressihle, ideal air- 
stream. The airfoils and their wakes are replaced by surfaces of discon- 
tinuity (vortex sheets), the interaction between the vortices being 
neglected. Each airfoil is performing both bending and torsional oscil- 
lations, t\TO adjacent airfoils being out of phase by a prescribed amoixnt. 

Under these conditions Euler's equations of motion are first linear- 
ized and then solved for the pressure distribution over the airfoils. 

The Biot-Savart theorem giving the induced velocity at any point on the 
airfoil due to the vortex field must also be used to obtain a solution. 
Once the pressure distribution is known the aerodynamic 'lift and moment 
acting on the airfoil can be obtained by integration. The con^ilete solu- 
tion is given in appendix B. The final equations and results will be 
given here. 


General case . - The aerodynamic lift and moment acting on an unstag- 
gered cascade as obtained in the manner outlined can be expressed in a 
form similar to that given in reference 10 for the isolated airfoil. 



where (all symbols are defined in appendix A) 


a elastic axis position as fraction of semichord b, measured from 
midchord, positive toward trailing edge 

b semichord 

h bending deflection, positive downward 

a torsional deflection about elastic axis, positive for increasing 
angle of attack 

p mass density of air 

(D circular flutter frequency 
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Ttie aerodynamic coefficients 1^, are functions of the 

following parameters: 

k reduced frequency, oSb/V 

Ertm phase angle between adjacent blades; m=0, l/2, l/4 
X cascade geometric parameter given by X = ut/s 


For the isolated airfoil these coefficients are functions of only the 
reduced frequency k- For convenience, the following additional param- 
eters are used: 


a = coth -g 
e = coth X 



where 

s spacing in units of semichord 
V free-stream velocity 

The following scheme, as outlined in appendix C, is now iised to calculate 
the aerodynamic coefficients: Let 

a© = 0; a^ = 1; bo = 1 
aj+l = i - 2^j] 3 

^j = '^j-1 ^j+1 



E 

i=0 


^i^i+J 

-2i 


'-j 


j 5^0 
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^j+1+21 °j-l-21 

^ i=0 (2i+l)a^^ 


= d 


= ^ *^j+l +21 '^j- 1-21 ^ 

1=0 ( 2 i+l)cj^^ 

m 

f ^ = CA + 2 2 _/ * 


'-J 


'i+1 


J >0 


J > 0 


(-1)'’ (c^.i + 2c j 


Dj = (-1)'^ (dj_i +2dj + dj,l)fj 
Ej ^ (-l)*^ (e^_^ + 2e^ + 

E, = (-l)*^ (d. , + 2d. + d . ) C‘^.1-21'*' ^.1-1-21^ 


Gj = (-l)*^ (d^^^ + 2dj + dj^3^) y ^ (®j_21 ®J-l-2i 


1=0 


Aq = 1 

= iS 


r(z) = gamma fimctlon of argument z 


(x,j) = x(x+l)(x+2) • * • (x+j-l) = 

F(i/ 2 , xj X + 1/2 + qj l/r)= ^ + 1 / 2 % q!i5(l! 


1=0 


±T 


Is the ordinary hypergeometrlc function. 
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% - E *; X + 1/2 + «i/t) 

g =0 


, , (x-m)^ 1/2-m r(x)r(l/2 1 ^ 

R 2 _ 1 + T l^ j Vp p ^ 

J=1 p=0 


J+p 


^ ,1/2- E E E ^ ®p 


i+j+p 


i=l j=l p=0 


/ ** ) 

^ (l ^0 + 


J =1 




’“' \ J=1 . 


P-(X,m) = - 


E 


8 e 

^3 4 
^ J =1 


46 

A. (J 


= + 


(l ®o + E 

\ J=i 


E. 


PgCX^m) = - ^ ^ V.J 

j=l 


C(k,X,m) = - ik + 


E 

d=i 

2ik + Kg 

~v 


P(k,\,m) = IkC - 


,2 21k - 2k + R„ 

k 3 
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The aerodynamic coefficients can now be obtained as follows: 



The geometric parameters X, e, cr, T are functions of the spac- 
ing and have been tabulated in table I for different spaeingjs. 

The quantities a^ and bj are functions of only the phasing m. 
The first 20 values of bj have been tabulated in table H for values 
of m equal to 0 and l/2. It is to be noted that for m = 0, bj is 
eq;aal to 1 for all j. 

In the equations for Cj, dj, and ej, the index J never taJtes on 
negative values. The quantities d_j, e_^ must always be obtained 

by calculating Cj, dj, and ej. 

The P-j^, Pg^ ^3^ P47 and P^ are functions of the cascade geometric 

parameter X and the phasing between blades m. The functions C and 
D are functions of X, m, and the reduced frequency k. It is shown in 
appendix B that the case m = l/4 (90° phasing) reduces to the case 
m = 1/2 (l80° phasing) with t-vrice the spacing. For the case of an iso- 
lated airfoil,- the functions reduce as follows: 

P3-.O 

P4 1/2 

P5 1/8 

C(X,m,k) -*■ C(k) Theodorsen's function 
D(X,m,k) - C(k) - ^ [1 - C(k^ + ^ 


s • 

X 0 


P3_-»- 2 
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The coefficients 1^, and then reduce identically to those 

of the isolated airfoil as given in reference 10. 

The values of .Aj are independent of geometry or phasing. The 
first 30 values are tabulated in table III. 


Specif case, m = 0 ; For the case of all the blades vibrating in 
phase tnt=0), a few sin 5 >lifications can be made. 

a^ = 0 


bj = 1 for all j 



^1- 


2 

He-1) 


Po = 


\2(e-i) 


la 


e ^-1 


„ .2-1 1 -Cv^V (1/2. x; *+1/2; 1/-^) 

C = K In 2 — + i72 

e T ' F(-l/2, Xj X + l/2; 1 /t) 

For this case, therefore, the bending coefficient Lj^ is obtained 

in closed form. The preceding results -can be obtained by summing the 
series for the appropriate quantities or by integrating 
ly of appendix B directly for m = 0. 

APPLICATION MB BESUITS 

As an example, the case in which all the blades are vibrating in 
phase (m=0) will be considered. Experimental data on compressors indi- 
cate that the only important type of vibration occurring is one of piure 
bending. For this case, then, the only aerodynamic coefficient that 
need be considered is 
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The functions P]_ and sxe plotted in figure 2 against the 

cascade geometry parameter X. The real and imaginary parts of the cir- 
culation function C are plotted in figure 3. It is to he noted that 
C is a function of hoth the spacing and the reduced frequency, -whereas 
the P functions are independent of the reduced frequency. 

In figure 3, over the range of reduced frequencies from zero to 0.6 
and spacings s = 1.67 to 3.33 (this corresponds to solidities of 0.83 
to 1.67), the real part of the circTilation function C is independent 
of reduced frequency for practical purposes, -whereas the imaginary part 
of C -varies linearly -with reduced frequency. 

With the P fimctions and the circiilation function C available, 
the aerodynamic bending coefficient \m.a confuted and plotted in 

figure 4. The -values for the isolated airfoil are also plotted on this 
figure. It is to be noted that (as well as the other coefficients) 
is complex. The real part of is in phase with the displacement and 

can do no work; the imaginary part is in phase -with the velocity and cor- 
responds to the danping component and can do work. Whether this system 
is stable therefore depends on the sign of the out-of -phase component or 
imaginary part of L^. The real part of is practically independent 

of reduced frequency and -varies only slightly -vrlth spacing in the range 
considered. The imaginary part of Lj^, however, varies appreciably 

both -with reduced frequency and spacing. This is further illiistrated 
in figure 5, where the coefficient has been plotted against spac- 

ing for a reduced frequency of 0.4. The results of figure 4 agree 
-vrell -vrlth those of reference 9 for a spacing s of 2.0, which is the 
only one considered there. Seme of the data from that reference are 
plotted in the figure. 

Since the imaginary part, or out-of -phase component of Lj^, corres- 
ponds to the aerodynamic damping, figure 4 shows clearly that the effect 
of cascading is to reduce greatly the aerodynamic danping when the blades 
are -vibrating in phase. For exanple, at a reduced frequency k of 0.4 
and spacing s of 2 (solidity of l), the aerodynamic damping is approxi- 
mately one -half the -value for the isolated airfoil. At a reduced fre- 
quency of 0.1, the damping is reduced to almost one-third the isolated 
airfoil value. The higher the solidity, the lower the aerodynamic damp- 
ing becomes . However, for this particular phasing (all the blades -vibrat 
ing in phase), the damping never actually goes to zero in pure bending 
except in -the limit of infinite solidity. - 


CCNCUJDING EEMARKS 

The equations for the oscillatory aerodynamic forces acting in sev- 
eral cases on a cascade of airfoils in potential flow are derived. The 
aerodynamic bending coefficient is calculated and plotted for the case 
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•where n.1 T the blades are ■vibrating in phase. For this case, the effect 
of cascading is to reduce greatly the aerodynamic dati 5 )ing force. 

The res'ults of this investigation can be used in the S'tudy of stall 
flutter phenomena in a cascade in a manner similar to the use of the 
classical flutter theory of the isolated airfoil for the stall flutter 
of the isolated airfoil. Aerodynamic time lags can be Introduced into 
or characteristic times (ref. 2) can be deduced from the aerodynamic lift 
and moment equations. The linear part of the blade characteristic (ref. 
U) is directly obtainable, and the nonlinear part, if tnown, might be 
treated by the methods of reference 2. The general usefulness of this 
approach, for both the isolated airfoil and a cascade, must, ho'wever, 
still be determined. 


Le-wis Flight Propulsion Laboratory 

National Ad'vi.soiy Committee for Aeronautics 
Cleveland, Ohio, July 16, 1954 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report: 

elastic axis position, measured from midchord, as a 
fraction of semichord b, positive toward trailing 
edge 




(k, X,m 
(k,X,m 


) 






f(t),f(ri),g(n) 

h 

1 

k 

L 




a 


M 

“a 


m. 


coefficients in recurrence formulas given In text 
semichord 

fxinctions defined in text 

coefficients used in evaluating integrals, given in text 

base of natural logarithms 

functions of indicated variables 

bending displacement of airfoil, positive downward 

integrals defined by eqs. (B40) and (B41) 

reduced frequency, cbb/V 
lift per TUiit span 

aerodynamic coefficients in lift equation, defined in 
text 

moment about elastic axis 

aerodynamic coefficients in moment equations, defined 
in text 

phase angle lag between any two adjacent blades, as 
fraction of 2it ladlans, 8/23t, appendix B 
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P’ 


s 

T 

t 

u 

V 

V 


X,Xj_ 

y 


a 

P 

r 

r 

5 

e 




functions defined in text 

pressure at any point on airfoil, function of time 
amplitude of pressure difference 

pressure difference at a point on airfoil, p^' - Pj’, 
fimction of time 

spacing tetween tlades in units of semichord h 

T-1 

geometric parameter defined by 

T 

time 

local velocity component in free-stream direction 
free-stream velocity 

local velocity con^ponent perpendicular to free-stream 
direction 

$ 

coordinate in free-stream direction in units of semichord 

coordinate perpendicular to free-stream direction in 
units of semichord 

angular displacement of airfoil 

stagger angle 

total circulation about airfoil 

vorticity distribution for reference airfoil 

vorticity distribution in wake of reference airfoil 

phase angle between two adjacent blades 

geometric parameter defined by c='coth X 

transformed variables of integration defined in appen- 
dixes C and D 


X 
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X 

P 

CT 

T 

<P 

OJ 


geometric parameter defined by — e 

s 

air density 

geometric parameter defined by a = cotb 

geometric parameter defined by T = 
velocity potential 
circular frequency of oscillation 



Subscripts : 

i,j,n,p,q summation indices 

u upper surface of airfoil 

2 lower surface of airfoil 

The subscript n is used also to indicate. the n^^ airfoil 




14 


MCA TN 3263 


C 


APPENDIX B 


LIFT MD MCMENT EQDAnONS 
General Theory 


The analysis of the oscillating airfoil of Infinite aspect ratio 
in a cascade proceeds along classical lines. The airfoil is assumed to 
he of small camber and thickness and performing infinitesimal oscilla- 
tions in an incon^jressible ideal fluid moving at a velocity V at 
infinity. 


The airfoil and its wake are replaced hy a surface of discontinuity 
or vortex sheet of strength r' over the airfoil and y^’ in the wake. 
This s'urface of discontinuity is assumed to lie in a horizontal plane 
parallel to the direction of flow. The vertical displacement due to the 
Interaction between vortices of this sheet is negjlected, this displace- 
ment being assumed sma.n coarpared with the horizontal motion. The 
asstaiption of small, perttirbations to the free-stream velocity permits 
the linearization of Euler's equations and, with. the introduction of a 
velocity potential, Bernoulli's equation for nonsteady Incompressible 
flow is obtained. The derivation 'is given here in detail for completeness. 

Euler's equations for two-dimensional flow can be written 




(V+u) + (V+u) (V+u) + V (V+u) 


1 dp' 


^ + (V+u) ^ + V = - i ^ 
^ dx ^ P oy 




(Bl) 


Considering only first-order terms and realizing V is a constant give 
the linearized Euler eq-uations as 


du , du 1 Sp' 

3t 

dv . dv 1 dp' 



The velocity potential (p 


is now introduced; 



(B2) 


(B3) 
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Eqiiation (B2) can now be written 




^ s (s) " ' p 


(B4) 




wbicb leads directly to Bernoulli’s equation for nonsteady motion 

a<p a<p , 1 , \ 


(B5) 


The airfoil and its 'v/ake are now replaced by a vortex sheet of strength 
T' on the airfoil and y ' in the wake. The difference between the 

W 

velocities of upper and lower surfaces of the vortex sheet is therefore 
Y*' or • For a point on the airfoil, therefore, 


& ~ ^ 


= yt 


■vdilch upon integration gives 


<p^ - (p^ = b 




(B6) 


Equations (B5) and (B6) then give 


( f X 

Si 


(B7) 


Continuity of pressure in the wake reqtiires that be zero there, and 

equation (B7) becomes, for the wake. 


, b a , 

^w + V ^ V 


. b ar> 

^1 - V ^ 


(B8) 


where 


■ 


dXn 


(B9) 
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The assiTTT^tlon of singale harmoixic motion 'will now he made, so that 


yi = ye^ 

V' = ve^ 




r' = r e 


tot 


T ' = T e 
' w ‘ w 


tot 


e 

Eqxiation (B8)' now becomes 


r • = re^ = e^ 

tot 


X dx 


J 


r + ih 

w 


r 


r„ dx3_ 


’ w 


ihr 


(BIO) 


(B8a) 


Eqmtion (B8a) is a simple nonhomogeneous integral equation with the 
kernel equal to 1. Its solution is 


= - itcre 


-ik(x-l) 


(BU) 


Equation (B11) gives the vorticity distribution in the wake as a ftmction 
of the total circulation around, the airfoil. It can also be obtained, 
from the condition that the total circulation around the system conipris- 
ing the airfoil and the wake must equal zero. 

With the use of eq-uation (B7) the pressure distribution over the 
airfoil can now be written in the following form: 




X 


r 


(B12) 


The lift and moment about the elastic axis can be foimd by inte- 
gfrating the pressure distribution as follows: 
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L = -pVb 


M = -pVb 


x: 





>(B13) 


Xi 


fO 

I 

g 


Vortlcity Distribution in a Cascade 

I Consider an infinite cascade of airfoils of chord Z, stagger angle 
P^-and spacing s as indicated in figure 1. ^Ehe airfoils are replaced 
by vortex sheets, the vorticity of the n'^^ sheet being designated by 
use of a subscript n. The Biot-Savart theorem is applied to give the 
induced velocity at a point x on the reference airfoil due to an element 
of vorticity of strength located at the point (x^, ns cos p) . 


dv_ = — 


-r^Cx - X;j^)dx3^ 


n. r 2 ^ 

2itj (x - Xj^) + (ns cos p) 

Measuring from the y’ axis, this can be written 

Tj^(x - Xj^ - ns sin p)dx^^ 


(B14) 


•dv^= - 


2«J(ns cos p) + (x - Xji^ - ns sin 


(B14a) 


Summing over all n and integrating frcm -1 to •• chaaage equation 
(B14a) to 




- a? 


L 


r^(x - - HB sin p)dx^ 


/-I 


^ l^ns cos p + i(x-X2-ns sin cos p - iCx-x^^-ns sin p)~ 


(B15) 


The assumption is now made that any two adjacent blades are out of phase 
by an angle 


6 = 2j0a 0 < m ^ 1 
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This angle is constant through the cascade. The vorticity distribution 
of the n^^ blade is therefore out of phase with that of the reference 
blade at the origin by an angle 

= 2ram 


or 


iBitmn 

r = r e 

n ' 


(B16) 


where y is the vorticity distribution for the reference blade. Siib- 

stituting r into equation (B15) results in 
n 


v(x) = - -57 



r(xi) 


se 




E 


iZrtmn 


ax, + 




E i 2 jtnm 


i(xn-x) 1 

n — " >■ " ' 

^ j a 


se 


ip 


se 


-ip 


n - 


se 


ip 


The infinite sums can be evaluated (see ref. 12).. 



-\(l-2m) (x-Xn ) 
e 

sinh X(x-X 2 ) 


(B18) 


X = i 
8 


Substituting now equation (B18) into equation (B17) gives 


vU) = 1 


r(xj^)e 


-X(2m.-l) (x-Xj^) 


sinh X(x-, -x) 




r(xi) 


E(2m-l)(x2-x) 


'-1 


/-I 


sinh X(x]^-x) 


dxi 


(B19.) 


- It . -IP 
X = - e 


It shoxild be noted that althou^ eqmtion (B18) is not valid for 
m = 0, equation (B19) is valid for all m. Equation (B19) gives the 
induced velocity at a point on the airfoil due to the ccm^lete vortex 
field of the cascade. It reduces to the well-known equation for t];ie 
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isolated airfoil wtien the spacing s goes to Infinity, If use is made 
of the condition that the normal induced velocity must equal the normal 
component of the velocity of the airfoil, then equation (B19) gives the 
relation between the normal con^ionent of motion of the airfoil and the 
voirticlty distribution. The problem then becomes one of solving equation 
(B19) for the vorticity distribution as a fxmction of the airfoil motion. 
Once the vorticity distribution is known, the pressure distribution can 
be calculated from equation (B12). 

The solution of the general equation (B19) will not be attempted 
here. Instead, several special cases of zero stagger will be considered. 
For this case, p = 0, X = X = it/s, and equation (B19) reduces to 


^ 



w 

r(xj)co6h X(l-2m)(x2^-x) 
sinh X(xj^-x) 


dxi 


(B20) 


Equation (B20) wi ll now be solved for the following three values of m: 


m = 0 
1 

m = 2 
1 

m = - 
For m = 0, 


all the blades are in phase 

-.o 


adjacent blades are 180 out of phase 
adjacent blades are 90° out of phase 


v(x) 


X 

2jt 


2^ 

2jt 



'-1 



cosh X(xn -x) 

r ( xi ) - -- ^7 r dx , 

sinh X^x^-x; 

• e^(xi-x) - 


(B21) 


The second equation follows from Kelvin’s theorem, since Y'(xj_)dxj=0. 

For m = 1 / 2 , \J-1 


v(x) 


2^ 

231 



(B22) 
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For m = 1 / 4 , 



cosh ^(xj-x) 
slTih X(x2^-x) 


dXn 


X/2 

2 jc 



rixj) 


sinh 


K^i-x) 

j 


dx. 


(B23) 


The case for m = l/4 therefore reduces to' the case of m = l/2 ■with 
t\7ice the spacing. The same result can he obtained from purely physical 
considerations. Thus only equations (B21) and (B22) have to he solved. 

In order to solve equations (B21) and (B22)^ a transformation of 
variables similar to that used in reference 6 'will he made. Let 


(B24) 


in terms of \i and tj, 


t(m) Ml r~ ^ P ajL 

- (^+>1) ^ (eSVXt-'ir'^ 


(B25) 



tanh Xx = p./e 

tanh Xxj^ = il/e 

■ tanh X = l/e ^ 

Eqmtions (B21) and (B22) can then he \nritten 
making use' of equation (B20), to-heccme 


•where m takes on the value of 0 or l/2. Or, more briefly. 

Si ^ 

■vdiere g(p.) and f(Ti) correspond to the appropriate parts of equation (B2S). 

Equation (B26) has an expJjLcit solution for f(Tj), first obtained 
in reference 13 and proven rigorously for the real domain in references 
14 and 15. With the condition f(l) finite, the solution is 
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The solution of equation (B25) can now he written directly. 


(B27) 


r(^) 

2 2 

t -T] 


(e- 


vin. 




m-1 


*1+4 (e-n) 
(e+4f 


n-4 


ikT it 
2jt ® 



(g-Tli) dT^i 




(B28) 


The second integral can he partly evaluated hy reversing the order of 
integration: 



tn-4)C4-ni) 


It 



(B29) 


Substituting equation (B29) into equation (B2B) results in 


m-1 


t(ti) =2.1:^ (c+Tl) 

62.^2 « 1+q 




m-1 


[1+tx (e-ii) 

(e-Hxr 


^ _ 

Tl-4 


itr ±k 
-TT- e 




(B30) 


also 

r = - r dq (B31) 

X J-l 

Equation (B30) can he integrated as Indicated in equation (B31) to obtain 
the total circulation around the airfoil f . 
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(B32) 


Equations (B30) and (B32) give explicitly the vorticity distribu- 
tion over the airfoil. With the vorticity distribution known, the pres- 
sure distribution and hence the lift and moment can now be ceiiculated. 


Lift and Moment on Airfoil in Cascade 


Consider the airfoil shown in figure 6. The upward displacement 
at any time t of a point x on the airfoil is 


y(x,t) = - h - b(x-a)a 


(B33) 


The induced velocity v(x) must satisfy the condition that the flow 
is everywhere tangential to the svirface. If the induced velocity in the 
x-direction is small compared with free-stream velocity V, this condi- 
tion leads to 


v(x) 




= - (h + Va - bad + box) 


or, in the t] coordinate. 


■v(ti) = -(h + Va - haa) - ^ a log 


e-n 


The lift and moment are given by 

a 


= b Zp(x)dx = ^{K)^K 




M = h 


s: 


Ap(x)(x-a)dx 




(B34) 


(B35) 


(B36) 
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Combining equations (B13), (B30), (B32), (B35)^ and (B36) and per- 
forming considerable algebraic manipulation, the lift and moment can be 
written in the following form: 


L = I + + a) l^Ja 

M = - I + I + [m^ - (I + + 1%) + (I 




(B37) 


tTbere 

I^ = PgCTt.m) -^P^(X,m) C(k,X,m) 




L^= I + PgCX^m) - i PgCX^m) 1 + C(X,m,k)J g” C(X,m,k) 

Lm 3c 


/ T A Pi (X,m) 

^ = ( I ■" SK ■ + —2 — D(X,m,k) 

\ ^ 3c 


‘a ^ i(l'a + \ - I ^ 


> (B38) 


M. 


PglX^m) 


|c(X,m,k) + D(X,m,k)J - ~ P^(X,m^ D(k,X,m) 
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P (X,m)= 


PpCXjia) = - 


2,2 2 
It X 


PgCx^m) = - 


P^^CX^m) = - 


P5(X,in) = - 


2it2x3 3 


4it2x3 


8ic2x^ 


C(X,m,k) = 


D(X,m_,k) = - 


e 2 

1 + Ik + k^ e L, 

2itX^ 


1 + ik - I k^ - k^ e^Q 

8icX 

^ ^ Xit ^ ® ^6 


and the I’s represent the following integrals: 

^ J-l J-l (£+>:)“ Vl-li n-li ^ 


^2^ 



di] 


(B39) 


t-\l Tl-Ji 


dTi 


Wi-ii n-ti 



dTi 


^ log liH ^ iJ 

1-p ^ 6-H 11-IJL ‘I 

(B40) 
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CO 

CD 

CD 


tJ) 

I 

& 

o 



I e| > 1 I III ^1 
Ini < 1 InJ > 1 

The integrals in equation (B40) are functions of only geometry and 
phasing between blades ( X and m) . ' The P functicxos are therefore 
functions of geometry and phasing only.- The integrals in equation (B4l) 
are functions also of the reduced frequency k. The C and D func- 
tions are therefbre functions of reduced frequency as well as of geome- 
try and phasing. The integrals in eq^tion (B40) are evaluated by inte- 
gration around a closed contour in the cQng)lex plane. The details are 
given in appendix C. The integral s in equation (B41) eire e-valuated by 
means of the hypergeometric integral. The deiails are given in 
appendix D. 


In the limiting case of the isolated airfoil (X = O), the following 
limiting values axe obtained for the various functlcms: 


\ 0 
C(X,m,k) 


Pl(X,m) 2 PgCx^m) ■* 1 

P4(X,m) -*■ 1/2 P5(x,m) 1/8 

C(k), Theodorsen’s function 
- C(k) - ^ (1 - C(k))- + ^ 


D(X,m,k) 



26 


MCA am 3263 


Equations (B37) and (B38) then reduce to those for the Isolated 
airfoil} and reduce identically to the coefficients 

for the Isolated airfoil given in reference 10. 
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APEEKDIX C 


EVALUATIQU OF INTEGRALS TO I5 

The Integrals I]_ to I5 are evaluated by integrating around a 

closed contour in the complex ^-plane, where the mapping function re- 
lating the t]-plane to the g -plane is given by 

■a=|(^ + iA) ■■■' ' ■ (Cl) 

By use of equation (Cl), equation (B40) can be written as double 
integrations around the closed contour C (fig. 7). Let 

= e + Ve^-l ^2 = ^ 

^3 = - e + ^4 = - e - = ~ ^1 


Then 


' ^ 4, (C-?i)"(C-Ca)" « 4 (C'-C3)=(C--C,r. ic'-55(c-Vc) . 

I 1 f r (f-q)°-^(C--Car^ ft.-uV 


2 (C'-CjUC'-Ci) 


* ^ 4 ' ,(?-cj)U-c«) 4 (c'-cjnc'-c^r 

^ ^ 4 u-^iAf-Cz)” ' L 4 (C'-Cjnc'-Ct)® (5'-0(t’-Vc) 


, (C'-Cj)(C‘-C.) 

W Tc,^c-)!c-- 4) 


is-t 


1 r (c-Ci)°-^(t-C4)"-^ «:-i)g r (ci-c)(c.cg)f r (f-Ci)°-^(c'-C3)°-^ (c.-m 

^4 ' L°® fC-C5)(C-C«)J J (C’-CjnC’-Ci)" (c'-0(c- 


(C'-C3)(C'-C«) 
VC) (Ci-C')CC'-Cg) 

(C2) 


It is to be noted that in integrating around the closed contour C, the 
integration path C* in the ii-plane is actually being traversed twice. 
Each integration around C must therefore be divided by 2. 

. If the integrations in (C2) are to be performed, the properties "of 
the various integrands must be considered. It will suffice to discuss 
one of them, since they are al 1 of the same genereil type. 
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Consider the integral ^31 defined by 


'■31 





log 


a'-^3)cc’-^) 

(C3) 


The integrand of has branch points at ^2 == ^3 “ “ 

and ^ ■well as simple poles at ^ and l/t- transfor- 

mation (Cl) maps eveiy point in the rj-plane in-to a pair of reciprocal 
points in "the ^ -plane. That part of C -which lies below the real axis 
is therefore the reciprocal of that part of C which lies above the 
real -axis, and every point which is outside the contour C has its 
reciprocal point inside the contoiar. The singularities of the integrand 
appear in reciprocal pairs, one of each pair lying outside the contour 
and one inside the contour. The poles at ^ and l/^ lie on the con- 
tour (see fig. 7). 


Considering the branch points, branch cuts are now made bet'ween ^2 
and ^2 and frcm and to infini-by as shovm. The integrand is 
now unique and single --valued over the path of integration. Furthermore, 
it can readily be sho-wn that since the integrand of I 21 is in-variant 
xonder the mapping ^ l/^ ' , (this evidently must be so because of (Cl) ) , 
the sum of the residues at ^ and l/^ must vanish identically. Con- 
sequently, any other contour such as equivalent to C but exclud- 

ing the points ^ and l/^, may be used. 


The Integrand of ^31 is therefore analytic in the annular region 
enclosed by the curves Cj^ and Cg as indicated in the figure. It can 
therefore be expanded into a Laurent series about -the origin -valid in 
this annular region. If is the coefficient of the l/^' term in 

this series, then ' 


Igj^ = 2jtlA_^ (C4) 

In order to expand the integrand into a Laurent series about the 
origin -valid in the indicated annul ur, . equations (C2) -will be -written 
as foUo-ws; 
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Each Integrand can now he obtained In the form of a Laurent series about the origin by expandlM each 
bracketed quantity (noting that each of these series will converge at least within the annulus) and 
multiplying the resixltant series together. In this way_, for exangile, the coefficient of 1/5 ' can 
be obtained and the inner Integration performed. This coefficient will, in general, be an Infinite 
series in The process then has to be repeated to carry out the second Integration. The final 

result will be a doubly or triply infinite series in c. The process Is evidently extremely long and 
laborious. However, it can be greatly singilifled by eng)loying a set of recurrence- formulas which can 
be obtained by inspection of the series, whereby each term of a given infinite series is obtained in 
terms of a previously given series. The calculation scheme is as follows: Let 


M 

to 
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8q = 0; = 1; = 1 

"j+i “ I - ^j] 


+ ^j+1 


Cj = a 




i=^) 


2i - --J 


( 21 + 1 ) 


oo 

E ^j+1+21 + ^J-1- 
/ ^ 21 


j-1-21 _ 

'J /-j , , 21 ^“>5 

i=0 (21+1) 


Cj = (-1)'^ (cj_^ + 2cj + Cj+i)^j ^ ‘^i+nj 
Dj = (-!)•’ (^j.l + 24j + * 2 ^ 

Ej = (-1)^ (ej_j^ + 2ej + + 2 ^ 

J-1 

Jj = (-1)^ (ij.i + 2ij + Vl> E t^J-2i ^J-l-21> 


1=0 

J-1 


Gj - (-1)'^ (ij_l + 2dj + ^ Cej_21 + ®j-l-21^ 


1=0 
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Then: 
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Let 


APPENDIX D 


EVADOAIION OF UnECSALS Ig TO Ig 






1/2 


6+T^. 


, 1/2 = ^ 

6-1 


T 


(Dl) 




X-V-, 


The order of integration is reversed and the logs are expanded in order 
that eqmtions (M) can he ■written as follows: 


1-x-m 


I,. = 


T X 


26 




- 1/2 


rV2(i_v)l/2(i_Tv)-‘“ 


'0 


>-l 


(l-T'v) dv dviL 


1-x-m _ 

^ T t: “ T 

■h ~ 2e. i 

.3=1 0 


ti / V 


- 1/2 r ,j-v2(i.vf/2 


(l-Tv)"’^(l-T'v)~^ dv dv^ + ^ Ig 


1-x-m 


= 


T X 


EE^ f 

i=l j=l ^0 ^0 


i+d- 1/2 


8 Ee 

(l-v)V2^_oi,)-ni(i_ni.v)-l dv dvj^ + I7 logT - f Ig 
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The inner integral is the same for all three, except for the expo- 
nent of V . This integral can he evaluated in terms of Appell*s hyper- 
geometric function in two variables of the first (refs. 16 and 17). 

r (l-v)^/^l-'Pv)"“(l-T*v)-^ d« = r(B-M/2)r(5/2) 

O r(n+2) 


F^(n+l/2;m,l;n+2;T,T') (P 3 ) 


r (n) Gamma fimction of argument n. 

The hypergeometric function F, is defined in terms of a double 
infinite sum in T and T'. 


m 

Fi(n+l/2jm,l;n+2;T,T*) = X/ ~ ^|n+ 2 '^|pT^^^ T^(T')^ (Bi) 


where 


(n,p) s n(n+l) (n+2) • • • (n+p-l) = 
(n,0) = 1; (l,n) = nl 


(D5) 


F^ can also he expressed an a single infinite sum of ordinary hyper 
geometric functions. 


m 


(D6) 


where F^-p,ljl4m^ 1 - is a polynomial given by 


q=0 

P 




q=0 


V ^fi - 

&7q) 1 r T y 

(D7) 
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Upon substituting back into equations (D2), the following integral re- 
mains to be evaluated: 

= 

“ ^ r C^;i/2 - f ’*! Ws-Hii i/x) 


(ref. 18) 

Equations (D2) can now be written as follows: 



(D9) 


where 
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r (i+j+p+2)r (1/2) 


with 


and 


Ao = 1 


= j±4±E+V2 


i+J+p +1 i+j+p +2 "^tj+p 


-p = -^ E c-)-^ p(x/.,x,xw^,iA) 

(DIO) 


q=0 


Let: 




R, 


= 1 + i2!.^r I ^/2-m TM Th 



4 B_ a? 


r(xti/2) 


J+p 


^ 


R, 


= l+iL^T^l/2-mrl2tM^ > B T^+J+P 

4 r tx+1/2) ^ ^ ^ IJ i+J+P P 


(DU) 


Then equations (B38) for C and D can he written as follows: 

, , . 2ik + Rp 

CCX,m,k) = - ik + g 


v2 2ik - 2k^ + R^ 
D = ik C - — £ 

2 Rt 


(D12) 
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TABLE I. - GEOMETRIC PARAMETERS 


[^=1 

■} t = coth. X, T 


T = CT 

= coth i] 

s 

X 

e 

CT 

T X 10^ 

T 

1.0 

3.1416 

1.0037 

1.0903 

293.27 

1.0000 

1.1 

2.8560 

1.0066 

1.1220 

92.434 

1.0000 

1.2 

2.6180 

1.0107 

1.1573 

35.312 

1.0000 

1.3 

2.4166 

1.0160 

1.1959 

15.876 

.9999 

1.4 

2.2440 

1.0227 

1.2372 

7.9398 

.9999 

1.5 

2.0944 

1.0308 

1.2809* 

4.3474 

.9998 

1.6 

1.9635 

1.0402 

1.3265 

2.5757 

.9996 

1.7 

1.8480 

1.0509 

1.3740 

1.6235 

.9994 

1.8 

1.7453 

1.0629 

1.4231 

1.0756 

.9991 

1.9 

1.6535 

1.0760 

1.4733 

. 74615 

.9987 

2.0 

1.5708 

1.0903 

1.5249 

.53585 

.9981 

2.1 

1.4960 

1.1057 

1.5774 

.39687 

.9975 

2.2 

1.4280 

1.1220 

1.6308 

.30253 

.9967 

2.3 

1.3659 

1.1393 

1.6850 

.23585 

.9958 

2.4 

1.3090 

1.1574 

1.7401 

.18787 

.9947 

2.5 

1.2566 

1.1763 

1.7957 

.15238 

.9934 

2.6 

1.2083 

1.1959 

1.8517 

.12565 

.9920 

2.7 

1.1636 

1.2162 

1.9085 

.10508 

.9905 

2.8 

1.1220 

1.2372 

1.9657 

.088957 

.9888 

2.9 

1.0833 

1.2588 

2.0235 

.076178 

.9869 

3.0 

1.0472 

1.2809 

2.0813 

.065934 

.9848 

3.2 

.9817 

1.3266 

2.1985 

.050747 

.9803 

3.4 

.9240 

1.3740 

2.3164 

.040292 

.9752 

3.6 

.8727 

1.4230 

2.4351 

.032811 

.9695 

3.8 

.8267 

1.4734 

2.5552 

.027299 

.9634 

4.0 

.7854 

1.5249 

2.6760 

.023138 

.9568 

5.0 

.6283 

1.7957 

3.2867 

.012345 

.9190 
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TABLE II. - t)j COEFFICIENTS 


i 





m=0 

II 

a 


m=0 

m=l/2 

0 

1 

1 

11 

1 

0 

1 

1 

0 

12 

1 

.2256 

2 

1 

.5 

13 

1 

0 

3 

1 

0 

14 

1 

.2095 

4 

1 

.375 

15 

1 

0 

5 

1 

0 

16 

1 

.1964 

6 

1 

.3125 

17 

1 

0 

7 

1 

0 

18 

1 

.1855 

8 

1 

.2734 

19 

1 

0 

9 

1 

0 

20 

1 

.1762 

10 

1 

.2461 





TABLE m. - Aj COEFFICIENTS 




J 


0 

1 

16 

0.008232 

1 

.2500 

17 

.007546 

2 

.1250 

18 

.006950 

3 

.07812 

19 

.006429 

4 

.05468 

20 

.005970 

5 

.04101 

21 

.005563 

6 

.03222 

22 

.005200 

7 

.02618 

23 

.004875 

8 

.02182 

24 

.004582 

9 

.01855 

25 

.004318 

10 

.01602 

26 

.004078 

11 

.01402 

27 

.003860 

12 

.01240 

28 

.003660 

13 

.01107 

29 

.003477 

14 

15 

.009963 

.009029 

30 

.003309 






Figsure 1. - Cascade geometry. 
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Figure 2. - Variation of and Pg with spacing. 
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Figure 3. - Variation of real and imaginary parts of C(X,k) with reduced 
frequency k for various values of spacing s. Phase lag m, 0. 











Aerodynamic bending coefficient 
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Figure 5. - Variation of bending coefficient Lj^ vitb spacing 
Reduced frequency, 0.4; phase lag m, 0. 


s. 
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Figure 7. - Contours in ^-plane. 
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